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Abstract 

This paper examines the well-posedness of the Stefan problem with a dynamic 
boundary condition. To show the existence of the weak solution, the original prob¬ 
lem is approximated by a limit of an equation and dynamic boundary condition of 
Cahn-Hilliard type. By using this Cahn-Hilliard approach, it becomes clear that 
the state of the mushy region of the Stefan problem is characterized by an asymp¬ 
totic limit of the fourth-order system, which has a double-well structure. This fact 
also raises the possibility of the numerical application of the Cahn-Hilliard system 
to the degenerate parabolic equation, of which the Stefan problem is one. 
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1 Introduction 

The Stefan problem is a well-known mathematical model that describes the solid-liquid 
phase transition. Among many results in the literature (for example, [171 ED ED ED and 
so on), the following enthalpy formulation of the Stefan problem with the Dirichlet-Robin 
type boundary condition was studied by the L^-framework in |15j : 

3u 

— -A(3{u) = g mQ:=(0,T)xf2, (1.1) 

where 0 < T < -boo is a finite time and C d = 2 or 3, is a bounded domain 
with smooth boundary T. The unknown u : Q —)■ M denotes the enthalpy and f3{u) is 
the temperature; : Q —)■ M is a given heat source. In the model of the Stefan problem 
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: M —)■ M is a piecewise linear function of the following form: 



r < 0, 

0 < r < L, 
L) r>L- 


( 1 . 2 ) 


ks, k 


i > Q represet the heat conductivities on the solid and liquid regions, respectively, 
and L > 0 is the latent heat constant. Consider the initial-boundary value problem for 
this kind of partial differential equation. The dynamic boundary condition is a sort of 
differential equation that represents the dynamics on the boundary T. As the condition 
for solving the partial differential equation, the time derivative is included, and this is 
well treated like, Dirichlet, Neumann, and Robin-type conditions for various problems. 
Under the dynamic boundary condition of the following form: 


d(3{u) 


-F d^f3{u) = gr on S := (0, T) x T, 


(1.3) 


the existence and uniqueness of fjl.ip was studied in nna using a subdifferential approach, 
where the symbol di^ denotes the normal derivative on T outward from hi; : S —)■ M is 
a given heat source on the boundary. For a more general setting, we can hnd the result 

in [3]. 

In this paper, the well-posedness of the Stefan problem with the following dynamic 
boundary condition is studied: 


du 

— + d^(3{u) - Ar(3{u) = gr on S, 
at 


(1.4) 


where the symbol Ap stands for the Laplace-Beltrami operator on F (see, e.g., [201 Chap¬ 
ter 3]). If we simultaneously consider fll.ip on the bulk hi and fll.3p or fll.4p on the 
boundary F, then the setting of fll.4p is more natural than that of fll.31) . 

The main idea of the existence result is to approximate the Stefan problem from 
the Cahn-Hilliard system. In 2009, Goldstein, Miranville, and Schimperna studied the 
following Cahn-Hilliard system: For e > 0 


du 

'm 


— A/i = 0 in Q, 


fi = —eAu + /3(m) -1- eTi{u) — f in Q, 
with a dynamic boundary condition of the following form: 

du 

— + dug. - App = 0 on S, 
dt 

g = eduU — eAru + (3r{u) + enriu) — fr on S. 


(1.5) 

( 1 . 6 ) 

(1.7) 

( 1 . 8 ) 


The unknowns m, /i : Q —>■ M represent the order parameter and chemical potential, 
respectively. Let us recall some basic concepts. The Cahn-Hilliard system is characterized 
by the nonlinear terms /9 -|- tt and /5p -|- vrp, which are some derivatives of functions W 
and lUp, respectively. Usually referred to as double-well potentials, for example, W{r) = 
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WY{r) = (1/4)(r^ — 1)^, in this case /5(r) = /9r(’") = and 7r(r) = 7rr(n) = —r for all 
r G M. For the details, see [g Eg. Now, taking /^r = /? as fll.2p and letting e ^ 0 
in fll.Gp . we obtain fll.ip as the limiting problem of fll.Sp with = P(u) — f and some 
modihcation of —A/ = g inQ and di^f — Arf = gr in T, with the setting that the trace f\^ 
of / coincides with /p. Then, from (II.8p . we can also characterize the dynamic bonndary 
condition fll.dp by fll.7p . 

The existence and nniqneness of the Cahn-Hilliard system fll.5p - fll.8p . as well as its 
asymptotic behavior, was treated by [19] in the case when e > 0. Following this, many 
related problems were treated in [2119]. Recently, the anthor extended the pioneering work 
of [ig to the more general case in which f3 and /dp are maximal monotone graphs. This 
included the singular case for the sub differential of the indicator function 

of the closed interval [—1,1], and when tt and vrp are Lipschitz continuous functions. 
The well-posedness of strong and weak solutions for the initial-boundary problem under 
appropriate assumptions was discussed in [TT] , where the essential idea comes from mm- 

The present paper proceeds as follows. In Section 2, the main theorem is stated. 
First, we prepare the notation used in this paper. Next, we recall a known result for some 
equations and dynamic boundary conditions of Cahn-Hilliard type (P)^ in Proposition 2.1. 
The main theorem is related to the convergence of Cahn-Hilliard systems (P)^ to the 
Stefan problem (P) with dynamic boundary conditions. In Section 3, we obtain a uniform 
estimate that is useful in proving the main theorem. To guarantee sufficient regularity of 
the unknowns for (P)^, we start from the approximate problem (P)£,a- After obtaining 
all necessary estimates, we correct similar uniform estimates for (P)^. A proof of the 
main theorem is given in Section 4. The strategy of the proof proceeds in a standard 
manner. Based on the uniform estimates, we consider the limiting procedure £ ^ 0. 
The main theorem is applied under a more general assumption for /3, namely that it is 
not only a non-decreasing piecewise linear function, as in fll.2p . but also some maximal 
monotone graph. Therefore, we apply the monotonicity argument. Finally, in Section 5, 
the uniqueness of (P) is proved by showing the continuous dependence for the given data. 

A detailed index of sections and subsections is as follows: 

1. Introduction 

2. Main results 

2.1. Notation 

2.2. Main theorem 

2.3. Approximate solutions to the Cahn-Hilliard system (P)^ 

3. Uniform estimates 

3.1 Uniform estimates for approximate solutions of (P)e,A 

3.2 Uniform estimates for approximate solutions of (P)^ 

4. Proof of the main theorem 

5. Continuous dependence 
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2 Main results 

In this section, the main theorem is stated. First, we recall a previous result from [11] 
that plays an important role in this paper. Next, we present the main theorem. 

2.1 Notation 

We use the spaces H := F^(r2), V := II^(Q), := F^(r), and W := H^{T) with the 

usual norms | • |h, | • |u, | • kr, I ' Iw ^’^cl inner products (lOu; (lOhf, (lOud 

respectively. Moreover, H := H x Hr, V := {{z,zr) E V x Vr '■ zr = z\^ a.e. on F}, 
and := H‘^{Vt) X Hereafter, we use a bold symbol z to denote the pair (^, ^r) 

corresponding to the letter. Then, H, V, and W are Hilbert spaces with the inner 
product 

{u, z)h '■= {u, z)h + {ur, zr)Hr all u,z E H, 

and the related norm is analogously dehned as one of V or W . Note that, if 2: G V, then 
zr is exactly the trace z\^ of 2 on F, whereas if 2 is just in H, then z E H and zr E Hr 
are independent. Dehne m : iT —)■ M by 

m{z) := —-— < / zdx + / 2 rdF I for all z G H, (2.1) 

l“l + TI IJn Jt J 

where |r2| := Idx and |F| := Jp IdF. The symbol V* denotes the dual space of V, and 
the pair {■,-)v*,v denotes the duality pairing between V* and V. Moreover, we dehne 
the bilinear form a(-, •) : W x V —)■ M by 

a{u, z) ;= / Vu ■ V zdx + / Vr^r ■ VrZrdF for all u, z E V, 

Jn Jr 

where Vr denotes the surface gradient on F (see, e.g., [201 Chapter 3]). We also introduce 
the subspace of H as Hq := {z E H : m[z) = 0} and Vq := V (1 Hq with their norms 
\z\ho •= \^\h for all 2 G Hq and \z\vo ■= 0 ( 2 , 2 )^/^ for all 2 G Vq- Then, we can dehne 
the duality mapping F : Vq ^ by 


{Fz, z)v*,Vo ■= a{z, z) for all 2 , 2 G Wq. (2.2) 

Using this, we can dehne the inner product in Vq by 

{z\,zi)y.:={zlF-'zt)v^y, for all zj, € VJ. (2.3) 

Then, we obtain the dense and compact embeddings Vq^ ^ Hq see m for 

details of this setting. These are essentially the same as in previous studies [101 ISSl [24] . 

2.2 Main theorem 

In this subsection, we dehne the weak solution for the Stefan problem with a dynamic 
boundary condition. Then, we give the main theorem. 
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First, we present the target problem (P) of this paper, namely the Stefan problem 
with the dynamic bonndary condition: 

du . r,, ^ „ 

— - A/f(M) = g a.e. m g, 

/3(tir) =/3(M)|r, -^ + d^(5{u) - = gr a.e. on S, (2.4) 

m( 0) = Mo a.e. in 12, Mr(0) = wor a.e. on F, 

where the prototype Stefan problem is formnlated by the setting fll.2p . In this paper, the 
target problem (P) will be formnlated in a more general setting; see Dehnition 2.1. 


Remark 1. The dynamic bonndary condition is arranged from three previons resnlts 
BEIE] regarding not only the Laplace-Beltrami operator, bnt also the time derivative. 
Actnally, they treated d(3{ur)/dt. In the case of fll.2p . the hrst condition fl2.4p implies 
that Mr is not necessarily eqnal to the trace u\^ of u. More precisely, we will obtain 
Mr = M|p except in the mnshy region {x G 12 : 0 < u{x) < L} of the Stefan problem. This 
is because we can only expect (3{u{t)) G V, but u{t) ^ V (weak solution). 

Throughout this paper, we assume that 

(Al) is a maximal monotone graph in M x M, and is a sub differential /3 = 5/3 of 
some proper, lower semicontinuous, and convex function /3 : M —)■ [0, +cxo] satisfying 
/3(0) = 0 with some effective domain D{P). This implies /3(0) = 0. Moreover, there 
exist two constants Ci, C 2 > 0 such that 

> Ci|rp — C 2 for all r G M; 

(A2) TT : h2(vr) = M ^ M is Lipschitz continuous. Moreover, 

|7r'(r)| < 1 for a.a. r G M; 

(A3) geL\0,T;Ho); 

(A4) uq := (Mo,Mor) G V with mg G intZ2(/3), and the compatibility conditions /3 (mo) G 
Li(12),^(Mor) e LHr) hold. 

For simplicity, we have assumed that the derivative of tt is bounded by 1 in fl2.6l) . It 
is sufficient to assume that vr G 1T^’°°(R). We now dehne a weak solution. 

Definition 2.1. The pair (u,^) of functions u G H^{0,T-,V*) nL^(0,T;FZ') and ^ G 
L^(0,T; W) is called the weak solution of (P) if 


(2.5) 

( 2 . 6 ) 


^G/3(m) a.e. in g, ^r e ^(^r), = ^|r 


a.e. on S, 
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and they satisfy 


(uft), z) + «(t), ZT)y, vr + / + f Vr^r(i) • Vr^rc^r 

’ ^ Jn Jr 

= / g{t)zdx+ / gr{t)zrdT for all z = {z,zy) G V 

Jn Jr 

for a.a. t E (0,T) with m(0) = uq a.e. in and Mr(0) = wor a-e. on F. 


(2.7) 


Our main theorem is now stated. 


Theorem 2.1. Assume (Al)-(A4) hold. Then, there exists at least one weak solution 
KO of{F). 

Moreover, we obtain the following continuous dependence for the given data: 


Theorem 2.2. Assume (Al)-(A4). Fori = 1,2, let ^ uj^ak solution of 

(P) corresponding to the data and Uq\ in which we assume that m(nQ^^) = 

Then, there exists a positive constant C that depends only on T such that 



for allt G [0,T]. Moreover, if fJ is Lipschitz continuous, then 

Jo 

< ^0(1 + T) + ^ (2.9) 

for all t G [0,T], where cy > 0 is the Lipschitz constant of {3. 

Note that the assumption = m{u^Q^) gives us that uFliyt) — u'^Atyt) G for 

all t G [0,T]. This theorem implies the uniqueness of the weak solution of (P). 


2.3 Approximate solutions to the Cahn—Hilliard system (P)e 

In this subsection, we state the approximate problem for (P). For this aim, we recall 
a previous result m for the equation and dynamic boundary condition of the Cahn- 
Hilliard type (P)^. This can be written as the following initial-boundary value problem 
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flTTOH - fimH : For e > 0, 

dti 

— — A/i = 0 a.e. in Q, (2-10) 

^ G /S(m), /W = —eAu + ^ + e7r{u) — f a.e. in Q, (2.11) 

du 

Mr = M|r, /^r = Ai|r, -^ + - Ar/Xr = 0 a.e. on S, (2.12) 

.^r G/^(Mr), jjir = — sA-pUr + ir + ^T^{ur) — fv a.e. on S, (2.13) 

m(0) = Uq a.e. in Mr(0) = Wor a.e. on F, (2-14) 


where / : Q —!■ M, /r : S —)■ M, mq : H ^ IR, and Mor : F —)■ M are given data. In particular, 
/ and /r are constructed by g and gp at a later point. In the main theorem, we treated 
the general setting (Al) and (A2) of f3 for the Stefan problem with some suitable vr. In the 
setting of fll.2p . one example of tt : M —)■ M is a piecewise linear function of the following 
form: 


/3(r) 


{ ksT if r < 0, 

0 if 0 < r < L, 

k^{r — L) if r > L, 


' L 
2 



if r < 0, 
if 0 < r < L, 
if r > L, 


(2.15) 


of course, in this case ^ 
(3 is obtained by 


(3{u) and .^r = I3{up). Therefore, 


(Al) and (A2) hold; actually 


^(r) 


hs 

—1 
2 

0 

ki, 


Lf 


if r < 0, 
if 0 < r < L, 
if r > L, 


so that (3 = d(3, and this quadratic function (3 satisfies fl2.5p . Because of this value of 
£7r, we can realize the double-well structure of the potential W = /3 -f- eif as the sum of 
primitives of (3 and stt. Therefore, (P)^ with the prototype setting fl2.15l) has the exact 
structure of the Cahn-Hilliard system. 


Remark 2. Viewed in terms of the Stefan problem, vr is an artihcial term. Therefore, 
we can prove the existence of the weak solution of the Stefan problem even if vr = 0. 
However, in this case, we no longer find the structure of the Cahn-Hilliard system (more 
precisely, the double-well structure) at the level of (P)e. If we construct the relationship 
between some equation of Cahn-Hilliard type and the original Stefan problem, we must 
choose a suitable vr depending on (3, which breaks the monotonicity of (3. This treatment 
is completely independent of the choice of boundary condition. We will focus on the 
convergence from the Cahn-Hilliard system to the nonlinear diffusion equation under the 
Neumann boundary condition in a forthcoming paper jl2j . 
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Here, we know that for each g := (( 7 , 5 'r) ^ Hq), there exists / G L^(0,T; Vq) 

such that 

/ '\/f-'Vzdx+ / Vr/r • Vr- 2 r<^r = / gzdx + grZrdT for all z G Vq, (2.16) 
Jo Jr Jo Jr 

a.e. in (0,T) (see fl2.23p . (12.241) . also [TT], Lemma C]). To state the previous known results, 
we use the change of variable v := u — mol, where tuq := m{uQ) and 1 := (1,1), 
the pair of constants. We also use vq ;= Uq — mol. We set /3(z) := (/3(z),/9(zr)), 
7r(z) := (7r(z), 7r(zr)) for all z e H. Under assumptions (Al)-(A4), [HI Theorem 2.2] 
leads to the following proposition of the existence and uniqueness for the equation and 
dynamic boundary condition of Cahn-Hilliard type (P)^ with respect to our f3 + en. 

Proposition 2.1. Under assumptions (Al)-(A4), there exists a triplet {v^, of 

V, G T; V*) n L°°(0, T; Wq) H L^{0, T; W), 

M,GL2(0,T;V), |,GL2(0,T;f/), 


that satisfy 


We{t),z)v,y^ + a{g,{t),z) =0 for all z e Vq, (2.17) 

= ea{v,{t),z) + (^^(t) +e7r(?;e(t) +mol) - f{t),z)^ 

for all z e V, (2.18) 

for a.a. t ^ with 

ie ^ ld{ve + mo) a.e. in 'Cr,£: £/d(ur.e + mo) a.e. on E. 

Moreover, '^ 5 ( 0 ) = Vq in Hq holds. 

We call the solution obtained by this proposition a kind of weak solution for the 
problem in fl2.10p - fl2.14p . This proposition is a direct consequence of [TT| Theorem 2.2]. 
Indeed, we assumed (Al)-(A4), and then, from the construction fl2.16p of /, we see that 
/ G L^(0,T; Vo). Thus, all of the conditions needed to apply [HI Theorem 2.2] can be 
corrected. 

To obtain uniform estimates with respect to £ G (0,1], we must consider the ap¬ 
proximate problem for (P)e, which is the same strategy used to prove Proposition 2 . 1 . 
Therefore, we give only a sketch of the proof here. For each A G (0,1], consider the 
approximate problem (P)£,a 

(''^£,a(^)> = 0 all zeVo, 

{t^e,xit),z)^ = \{v'^^^{t), z)^ + ea{ve,x{t), z) 

+ +mol) +e7T[ve,xit) + mol) - f{t),z)^ for all z G V, 

(2.19) 
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with i’£,a(0) = Vq in Hq, where (3\{z) := (3x{zr)) for all 2 : G iT and (3x is the 

Yosida approximation (see, e.g., [Ill5l|22]). Namely, /9 a : M ^ M along with the associated 
resolvent operator Ja : M —)■ M are dehned by 

Px{r) := “ '^a(?")) ■= — {I + •^/^)~^(^)) fo^' fol r G M. 

A A 

Moreover, the related Moreau-Yosida regularization /9 a of /9 : M ^ M fulhlls 

Mr) ■= taf - y + ?(o} = - •^A(r) t + ?(^A(r)) = (2.20) 

for all r G M. We know the basic property 

0 < (3x{r) < (3{r) for all r G M. (2-21) 

The problem (P)e,A can be solved (see, e.g., [lOl [HI [IH |23l |2l] ); more precisely, there 
exist v^^x £ hr^(0, T ; i?o)nL°°(0, T; Wq) nL^(0, T; W) and n^ x ^ that satisfy 

+ sdip{v,^xit)) 

= P{-I3x{ve,x{t) + mol) - en(ve,x{t) + mol) + f{t)) in Hq (2.22) 

for a.a. t G (0,T) with 'y£,A(0) = Vq in Hq, where </? : Ho —)■ [0, +cxo] is a proper, lower 
semicontinuous, and convex functional 

[ \Vz\^dx + l- [\VrZr\^dT if 2 G Wq, 
ip[z) := Jn 2 ip ('2.23) 

I+CX 0 otherwise. 

Here, the subdifferential dcp on Hq is characterized by dip{z) = {—Az, di,z — Ap^p) with 
2 G D{dip) = WnVo (see, e.g., HD Lemma C]). We also note that 

299 ( 2 ) = a( 2 , 2 ) = {Fz, z)v*,Vo = \^\vo z eVo- (2.24) 

Moreover, P : H ^ Ho is a projection defined by Pz := 2 — m(2)l for all 2 G H, and 
it satishes 

(20, Pz)ho = (zq, z)h for all Zo G Ho and z e H, 

|P2|vo = I^Iko a \z\v for all 2 G V. 

The standard strategy obtains a priori estimates with respect to A G (0,1] and considers 
the limiting procedure A —?■ 0. 

Remark 3. Let := + mol. Then, HD Remark 2] means that fl2.17p implies 

{u'^{t),z) +{u'^^^{t),zr)y,y + / VHeit) ■ Vzdx + / Vr/ar,£(t) ■ Vr^pdL = 0 

for all 2 G W. Moreover, thanks to the regularity G L^(0,T; VF), we see that (I2.18p 
implies 

C£ G /3{ue), fie = -sAue + ^£ + e7i{u^) - f a.e. in Q, 

^r,e ^/3{ur,6), fir,s = £di,Ue - eArUr,e + ^r,6 + £'^r{ur,e) - fr a.e. on S. 
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3 Uniform estimates 

In this section, we obtain the uniform estimates needed to prove the main theorem. 


3.1 Uniform estimates for approximate solutions of (P)e,A 

To prove the main theorem, we will use the uniform estimates that are independent of e 
for the solutions of (P)^. For this, we must start from the approximate problem (P)e,A- 
In this subsection, we obtain uniform estimates for the approximate solutions of (P)e,A- 
We recall the change of variable a := lie ,a + mol. 

Lemma 3.1. There exist positive constants Mi and M 2 , independent of e & (0,1/4] and 
X G (0,1], such that 

Al-Ue.AWl^g + < Ml, 

\v,^x{s)\y/s + 2 |^A(w,A(s))|^i(^)ds + 2^ |^A(Mr,e,A(s))|ii(r)C?S < M 2 

for all t G [0, T]. 


Proof. We test 02.221) at time s G (0, T) by ti£,A(s) G Uq which is considered in problem 
(P)£^a- Then, using 02.3p . we have 

A«a(s), W,A(s))jf^ + «a(s), W,a(s))^* + £(a(p(t;£,A(s)), W,a(s)) 

+ (1^/3aK,a(s)+ mol),'ii£,A(s))jj^ = {f{s) - eP7z{v,^xis)),v^^xis)) (3.1) 

for a.a. s G (0,T). Now, from the definition of the subdifferential, we have 


{PPx (w,a(s) + mol), W,a(s)) 

= / (3x{ue,x{s)){ue,x{s) - mo)dx + / /3A(Mr,£,A(s))(Mr,e,A(s) -mo)dr 
Jn Jr 

> [ P\{u^^x{s))dx - j3x{mo)\U\ + f ,5x(ur,f,A(s))dr - ,5x('mo)|r 


(3.2) 


for a.a. s G (0,T). Next, recalling the fundamental property of the chain rule with 02.61) 
and using the Young inequality, we have 


e{vs,x{s),PTT{ve,x{s) +mol)) 


V*o,Vo 


< -\v^^x{s)\y, + / \V 7 r{v,^x{s) + mo)\^dx + / |Vrvr(nr,£,A(s) + "lo) l^c^P 




(3,3) 


W,a(s),/(s)) 


1 

W,Yo - 4 



2 

Vo 


and 


(3.4) 
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for a.a. s G (0,T). Thus, correcting fl3.ip - fl3.4p . recalling the dehnition of the subdiffer¬ 
ential and using fl2.2ip . fl2.24p . and </ 3 ( 0 ) = 0, we see that 


+ 2|/5a(m£,a(®)) |l1(q) + 2|/3A(Mr,£:,A(^)) iLi(r) 
< 2 (|fl| + |r|)^(mo) + |^^£,A(5)rv5 + 2|/(s)rvo 


2 

Ho 


(3.5) 


for a.a. s G (0,T). If we take e G (0,1/4], then e — 2e^ > e/2 > 0. By virtue of the 
Gronwall inequality, we obtain 


- (I'^oIho + I'^olvs + 2T(|n| |r|)/3(mo) + 21 / 1 ^ 2 ( 0 ,t;Ko)) =• 


for all t G [0,T]. Next, integrating fl3.5p over (0,t) with respect to s, we obtain 

\veAs)\yo^s + 2 \A{ue,x{s))\^,^^^ds + 2 |3A(Mr,s,A(s))|^i(P)ds 

< I'^ol^o + l^^olvs + 2T(|G|-t-|r|)/ 3 (mo)2|/|^2(o^r;Ko) + [ |'^eA('S)|v£*^'® 

Jo “ 

< Mi(l + T) :=M2 

for all t G [0,T]. Thus, we have the conclusion. 


□ 


Lemma 3.2. There exists a positive constant M 3 , independent of e & (0,1/4] and A G 
( 0 , 1 ], such that 


2A 



V 


e,xiAv*/s + e\v,At)\ 


2 

Ho 




(3,6) 


for all t G [0, T]. 



(3.7) 


Proof. We test (12.221) at time s G (0, T) by v'^ ^{s) G Hq. Then, using the same method 
as for (13.3p . we have 


-^« a ( s )3 P'^AeA^) + )v*,Vo 



2 

Ho’ 


«a(s),/(s)> 


1 

HS,Ho - 4 


V 


£,X V 


(^)rvs+i/(^)i 


2 

Ho 


and 
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for a.a. s G (0,T). Then, using fl2.3p and the chain rule, we deduce 




d f d f 

— / P\{vs,x{s) + mo)dx + — / ^A(Tr,£,A(s) + mo)dr 
«/ w r 




I Vo 


(3.8) 


for a.a. s G (0,T). Integrating fl3.8l) over (0,t) with respect to s, we can hnd a positive 
constant M 3 , depending only on |no|vo, |^(mo)|l 1 (u)) l^(wor)|Li(r), M 2 , and |/|l 2 ( o r;Vo) , 
such that the aforementioned estimate fl3.6p holds. Next, to obtain fl3.7p . we recall (12.191) . 
(see also Remark 3). We have 


Me,A(s) = + A(“t,A(s)) +£■"•( 1 * 1 ,a(s)) - f(s) in V, 


(3.9) 


for a.a. s G (0,T). Therefore, the evolution equation (12.221) is equivalent to 


+ F(P/x^;^(s)) = 0 in V*Q, for a.a. s G (0,T), (3.10) 

with (13.9p . We test (I3.10p by G Vq, and integrate the resultant over ( 0 ,t) with 

respect to s. Then, using (12.21) and the Young inequality, we have 

^ Ws,xis)\v,ds + ^ 

for all t G [0,T]. Thus, using (I3.6p . we obtain (13.7p . □ 


Lemma 3.3. There exist two positive constants M 4 and M 5 , independent of e G (0,1/4] 
and A G (0,1], such that 

\Ue,x{t)\]j < M 4 , < M 4 , (3.11) 

|/?a(w,a(' 5)) |/?A(Mr,£,A(-s)) ^ ^5 (3.12) 

for all t G [0, T]. 


Proof. From (12.5p in (Al) and (I2.20p . 

^A(r) = ^|r - JA(r)f +^(jA(r)) 

> — |r — Jx{'f')\‘^ + Ci| JA(r)|^ — C 2 for all r G M. 

2a 

Therefore, from (13.61) in Lemma 3.2, we obtain 

^|w,A(t) - Jx{Ue,x{t))\'^H + ^ + C2|fi|, 

^|Mr,£,A(^) - </A(Mr,£,A(t)) 1^^ + Ci| JA(Mr,£,A(t)) < ^ + calT 
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for all t G [0,T]. Thanks to this A G (0,1] implies that there exists M 4 > 0, depending 
only on Ci, C2, M 3 , |r 2 |, and |r|, such that 

< ‘^\Ue,xit) - Jx{Ue,x{t)) + 2\jx{u,^x{t)) 

+ 2|Mr,£,A(^) - -^A(Mr,£,AW) l^r + ^| Ja (Mr,£,A W) l^r 

< M 4 

for all t G [0,T]. Then, there exists M 4 > 0, depending only on M 4 , |mo|, |fl|, and |r|, 
such that 


\ve,xit)\H ^ 2|n,,A(t)|^ + 2|mo|^(|fl| + |r|) 

< M 4 


for all t G [0,T]. Thus, (Id.llh has been obtained. Next, recall the useful inequality [IHl 
Section 5]. Indeed, from assumption (A4), we have rriQ G intT)(/3) (which is the criterion 
needed to apply the inequality), meaning there exist two constants C 3 , C 4 > 0 such that 


(3x{r){r - mo) > Csf3x{r) - C 4 
for all r G M and A G (0,1]. Then, 


(3.13) 

can be improved by 

{Pf^x K,a(s) + mol), ^^£,a(5)) 

= / Px{ue,x{s)){us,x{s) - mo)dx + / l3x{ur,e,x{s)){uT,e,x{s) -mo)dr 

Jn Jr 

>C3 [ ( 3 x{Ue,x{s))dx-C 4 \n\+C 3 f /^A (Mr,£,A('5)) - C4 |r| 

Jn Jv 

for a.a. s G (0,T). Therefore, by using fl2.3l) and the monotonicity of dip, 

C3 / IJx{Uep{s))dx + C 3 / /?A(Mr,e,A(^))^r 

Jn Jv 

< C4(|fi| + |r|) + A|n'_A(s)lHol^"'^(®)|//o + 

+ s\7T{u,,x{s))\jj\v,,xis)\jj^ + \fis)\^Jv,,x{s)\H^ 
for a.a. s G (0,T). Now, squaring the above and using Lemma 3.1 and fl3.1ip . we obtain 



|2 


ycs J /3x{ue,x{s))dx + C 3 

< 5cl{\Q\ + |r|)^ + + 5M4|<,(s)1 

+ 5£^M4|7r(we,A(s)) 1^ + 5M4 |/(s)| 


2 


Ik 


(3.14) 
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for a.a. s G (0,T). Thus, from (12.61) and (13.lip , we have 
|7r(we:,A('5)) |L = K(^£,a(' 5)) “ 7r(mol) + 7r(mol)|^ 

<2|y (|t£,a(s)|^ + |7r(mo)f)da; + J (^|ur,£,A('5) T + k("^o) 

< 2 ||ue,A(s)l^^ + | 7 r(mo)f (|f2| + |r|)| 

<2|M4+|7r(mo)f(|f2| + |r|)} =: Mg (3.15) 

for a.a. s G (0,T). We integrate the resultant (I3.14p over (0,t) with respect to s. From 
Lemma 3.2, there exists a positive constant M 5 , depending only on C 3 , C 4 , T, |f2|, |r|. Mi, 
M3, Ml, Ms, and \f\L'^{o,T;Ho) and independent of e G (0,1/4] and A G (0,1], such that 
dim holds. ’ ’ □ 

Now, recalling 03.Op and using the fact that dip{v^^xit)) G Hq for a.a. t G (0, T), 

we obtain 

+ £7^{ue,x{t)) - fit)) for a.a. t G (0,T). (3.16) 


Lemma 3.4. There exist two positive constants Mg and Mj, independent of e & (0,1/4] 
and X G (0,1], such that 

< ^6, 

M£,a(s)|^c(s < Mr 

for all t G [0, T]. 




Proof. From (12.ip and (I3.16p . we have 

h(/x,^^(s)) I' < { |/^aKa(s)) + £2|7r(n,,A(s))+ |f2| |/(s)|^ 

+ \/3x{ur,e,x ('®)) Lpr) + I T(nr,£,A (®))lLi(r) + |r||/r(s)|^^} 

for a.a. s G (0,T). Then, by integrating over (0,T), it follows that there is a positive 
constant Mg, depending only on T, |f2|, |F|, M4, Mg, |7r(mo)|, and \f\L^{o,T;Ho), such that 
the hrst estimate holds. Next, from the Poincare inequality (see, e.g., [m Lemma A]), 
we know that there exists a positive constant Cp such that 

I |2 ^ I |2 

Cppiv < l^lvo 


for all z G Mo- 


(3.17) 
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Therefore, 



^^e,\is)\vds < 2 



P ^^ e ,\ i ^)\ v^S + 2 





2 

V 


ds 


+ + |r|) ^ \m{fi^^^{s))\‘^ds 

<‘^ + 2{\n\ + \T\)Me =: M 7 

Cp 


for all t E [0,T], 


□ 


Lemma 3.5. There exists a positive constant Mg, independent of e E (0,1/4] and A G 
(0,1], such that 


j |/3A(“e,A(s))|H<*S < 


(3.18) 


for all t E [0, Tj. 


Proof. We test fl3.9p by f3xius,x{s)) E V, then 

e{d^{ve,x{s)), f3x{ue,x{s))) ^ + \(3x{ue,x{s))\]j 
< - >^v'^,x{s) - S7r{u,^x{s)) + f{s), f3Aue,xis))) 

for a.a. s E (0,T). Here, from the monotonicity of /3x, we have 
£('99?(w,a(s)),/3aK,a(s)))jj 

= e {-Ave,x{s))f^x{ue,x{s))dx + s / di,Vs,x{s)^^\{ur,e,xis))dT 

Jn Jr 

+ e j {-ArVr,eAs))l3x{ur^e,xis))dT 

= e j /3^Me,A(s))|Vn£,A(s)rc?a: + j /3^Mr,£,A(s)) | Vrnr,£,A( 5 ) l^c^r 
> 0 


for a.a. s E (0,T), where we have used the fact that /9A(M£,A(®))|r ~ fJx{ur,e,x{.^)) 
r. Moreover, by the Young inequality. 


(a^£,a(^) “ + fis),f3xAeAs)))jj 



+ 2A^|<a(s)|L + 2£2|7r(w,,A(s ))\h + 2|/(s)| 


2 

Ho 


for a.a. s E (0,T). Let us combineing the above inequalities, and integrateing in over 
(0,T) with respect to s. Then, recalling (I2.2ip and (I3.15p . and using Lemmas 3.2 and 
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3.4, we obtain 



(3^{u^^x{s))\]jds <4 + 4:\^ K ;,(s)|^^c?s 

+ 4 ^^ f \7z{u^^x(,s))\]jds + 4 [ \f{s)\lj ds 

Jo Jo 

< 4Mj + 4M3 + 4TM5 + 4|/|^2(o,T;iLo) 


for all t E [0, T]. 


□ 


3.2 Uniform estimates for approximate solutions of (P)e 

In this subsection, based on the result of Proposition 2.1, we obtain uniform estimates for 
the solutions of (P)e. Actually, by obtaining additional uniform estimates (see [11]) and 
passing to the limit in the approximate problem (P)e,A as A —)■ 0, Proposition 2.1 can be 
proved. Thus, for each e G (0, 1/4], there exists a triplet of 

V, E H\0, T- V*) n L°°(0, T; Uq) H ^^(o, T; tf"), 

/x,gU( 0,T;U), CGL^(0,T;ff), C e/3(n,) inL2(o,T;if) 

that satisfy fl2.17p . fl2.18p . and l’£(0) = Vq in Hq. Here, we also have 

u, E N^(0, T; V*) n L°°(0, T; V) n T; W), 


with the relation + rngl. On taking the limit A —>■ 0, we obtain the same kind of 

uniform estimates as obtained in the previous lemmas, namely 


[ \v'e{s)\v*ds < M 3 , (3.19) 

Jo 

e\v,(t)\l^<M„ (3.20) 

\u,(t)\%<M,, (3.21) 

< Mt, (3.22) 

[ \^i,(s)\lds<M^, (3.23) 

Jo 

[ |4(5)|^ds<M8 (3.24) 

Jo 


for all t E [0,T]. Using these estimates, we can prove the main theorem. 


4 Proof of the main theorem 

In this subsection, we conclude the proof of the main theorem by passing to the limit in 
the approximate problem (P)^ as e —)■ 0. 
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Proof of Theorem 2.1. By using the estimates fl3.19p - fl3.24p stated in the previous 
section, there exist a subsequence with —)■ 0 as /c —)■ +cxo and some limit func¬ 

tions ir G H\0,T;V*)nL°^{0,T; H), v G H\0,T;Vl)nL°°{0,T; Hq), G L‘^{0,T;V), 
and ^ G L^(0,T; H) such that 


—)■ V weakly star in L°°(0,T; Hq), 

(4.25) 

xxgj. —7- tx = n -|- mol weakly star in L°°(0, T; iT), 

(4.26) 

in L°°(0,T; Vo), 

(4.27) 

—)■ v' weakly in L^(0,T; Vq), 

(4.28) 

—)■ u' weakly in L^(0,T; V*), 

(4.29) 

M weakly in L^(0,T; V), 

(4.30) 

weakly in L2(0,T;ff), 

(4.31) 

£fc7r(xx£j-G 0 in L°°(0,T; iT) 

(4.32) 


as /c —)■ -l-cxo. From fl4.25p . fl4.26p . fl4.28p . and fl4.29p . the well-known Ascoli-Arzela theo¬ 
rem (see, e.g., [26l Section 8, Corollary 4]) gives 

mC(|0,T];W), (4.33) 

mC(|0,r];W) (4.34) 

as fc —)■ -|-cx). Now, integrating fl2.18p over (0,T), we hnd 

[ ^dt 

Jo Jo Jo 

+ ek [ {n{ue^{t)),'n{t))^dt- [ {f{t),'n{t))^dt 
Jo Jo 

for all r/G L2(0,T; V). (4.35) 

Letting /c —)■ cxo, we can use fl4.30p . fl4.27p . fl4.3ip . and fl4.32p to obtain 

[ for all r/G L^(0,T; V), 

Jo Jo 

namely, /x = ^ — / in L^(0,T;i/). Additionally, the information /x-|- / G L^(0,T; V) 
gives the regularity of ^ G L^(0,T; V), namely .^r = ■Clr Next, we take r] := 

xxe^ G L^(0,T; V) in fl4.35p . Then, using the positivity 


Sk a{v^i^{t),u^^{t))dt = Sk 


\v 




{t)\ydt > 0 , 


and recalling fl4.26p . fl4.30p . fl4.32p . and fl4.34p . we have 


limsup / {^^^{t),u,^{t))^dt < 
k —^-|-cx) «y 0 


it),^J'{t))y,ydt+ / {f{t),u{t))^dt 

uv ’Jo 
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Thus, applying [U Proposition 2.2, p. 38] with fl4.25p and fl4.3ip . we can deduce that 

^ G fi{u) a.e. in Q 

from the maximal monotonicity of jS. Namely, we obtain ^ G f3{u) a.e. in hi, ,^r G P{ur) 
a.e. on P. Finally, integrating fl2.17p over (0,T) with respect to t and letting k —)■ +cxo, 
we have 



Next, using (12.2p . we deduce 

v\t) + F{Pfi{t)) =0 in Vo, for a.a. t G (0,T), 

namely 

(n'(f), z) = a{f{t), z) for all 2 G Vo- (4.36) 

Therefore, recalling (12.161) and [HI Remark 2], we hnally obtain 


+ {u'j^{t),zr)y, y^ + / V^(t) ■ Vzdx + / Vr,^r(t) ■ Vr^rdT 


'Vr*,Ur 


= / g{t)zdx+ / gr{t)zrdT for all 2 = (2, 2r) G V, 


Jn Jr 

for a.a. t G (0,T), with m(0) = uq a.e. in hi and Mr(0) = uor a-e. on P. Thus, it turns out 
that the pair (u,^) is a weak solution of (P). 


5 Continuous dependence 

In this section, we prove the continuous dependence of the data. This theorem also 
guarantees the uniqueness of the component u in the solution. 

Proof of Theorem 2.2. For i = 1, 2, let pg ^ weak solution of (P) correspond¬ 
ing to the data Set ttiq := 771 ( 1 ^ 0 ^^) = From the weak formulation 

fl4.36p of the Stefan problem (P) for nb) = — ttiqI, we obtain 

+ for all 2 G Vo, (5.37) 

for a.a. s G (0,T). Now, = ^‘'*^(s) — m(£*(s))l G Vq. Therefore, we see that 

a{^^^\s),F-^Zo) =a{F-^Zo,P^^^\s)) 

= (2o,PCH^)>VS,Vo 

= (2o,PCn«))Ho 

= (2o,CHs))^ for all 2o G ifo C V*. 


(5.38) 
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Thus, we take the difference between equation fl5.37l) when i = 1 and when i = 2. Then, 
we set 2 := E Vq, and use fl2.2p to obtain 

< (nW(s) - n(2)(s), /(^)(s) - 

< - f^^\^)\vo 

for a.a. s G (0,T). Now, in the variational formulation fl2.16p used to construct /, taking 
2 ;= /, we obtain 



a.e. in (0,T). Therefore, by using the Young and Poincare fl3.17p inequalities, 

< 5|/“’(o) - + i {!/“'(<!) - 

< - i|/"w - 

£ 5|||9''’(s)-9‘"’(s)lk (5.40) 

for a.a. s G (0, T). Then, from the monotonicity of /5, we can apply the Gronwall inequality 
to deduce fl2.8p for some positive constant C that depends only on Cp and T. Here, we 
must take care to ensure that = wb) _ -u(2)_ If is a Lipschitz continuous 

function with Lipschitz constant C /3 > 0, then ^ = (/3 (m),/ 3(Mr)) and, with the help of the 
monotonicity of /3, we have 

< C;3(w(^^(s) (5.41) 

Thus, combining fl5.39p . fl5.40p . fl5.4ip . and fl2.8p . and integrating the resultant quantity 
over (0,t) with respect to s, we obtain 

° Jo ° Jo 

< C{l + T) |c/W(s) 


for all t G [0,T]. That is, we have obtained (12.9p . 


□ 
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